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Abstract 

We introduce a new tool for "unmasking" the composition of neutron star (NS) interiors which 
is based on the fact that the state of matter at high densities determines the statistics of both 
NS observables, the temperature-age (TA) data as well as the mass distribution. We use modern 
cooling simulations to extract distributions of NS masses required to reproduce those of the yet 
sparse data in the TA plane. By comparing the results with a mass distribution for young, nearby 
NSs from population synthesis we can sharpen two NS cooling problems. The direct Urea (DU) 
problem consists in a narrowing of the NS population at the mass value for which the DU process 
as the most effective cooling mechanism in the hadronic layer of the star can occur. The Vela mass 
problem is a broadening of the population beyond the range of the typical mass window of 1.1 - 
1.5 Mq. Applying this tool to modern EoS we discuss examples for pure hadronic stars which are 
in conflict with these constraints while hybrid stars with a color superconducting quark matter 
core can predict a satisfactory mass distribution, provided the smallest diquark pairing gap has a 
properly defined density dependence. 



1 Introduction 



The quest for understanding the state of nuclear matter at high densities beyond saturation at n > 0.16 
fm -3 has triggered a multitude of theoretical and experimental investigations. It is of principal interest 
to clarify the mechanism how the quark and gluon substructure of nucleons might manifest itself under 
extreme compression pp. One of the central questions is that of deconfinement: nature might choose to 
excite more massive degrees of freedom of the hadron spectrum such as hyperons and nucleon resonances 
in order to satisfy the Pauli principle before the quark-gluon plasma shall emerge as expected from the 
asymptotic freedom property of QCD at ultimate densities. Due to the presence of strong attractive 
correlations in this intermediate density regime, the transition of hadrons from well-localized quark 
bound states to delocalized, unbound states is expected to pass a transitory phase of strongly correlated 
few-quark complexes with frequent resonant scattering processes [21 [3] . Thereby a very rich structure of 
quark condensates can occur characterizing color superconducting phases in the low-temperature QCD 
phase diagram, most relevant for compact star physics [HIUE]. 
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Besides terrestrial experiments with relativistic heavy-ion collisions the interior of neutron stars 
provides a "laboratory" where matter under conditions of extreme densities occurs [5J. The problem is 
to identify the composition of neutron star interiors from their observable properties like, e.g., masses, 
radii, rotational and cooling evolution |7, ;"$]• 

Until recently the compactness of objects has been discussed as a characteristic feature of stars with 
a quark matter interior (strange stars) when radii are less than ~ 10 km [9J [10]. The situation has 
dramatically changed with the advent of precision measurements of high masses for objects like the 
pulsar PSR J0751+1807 [H] with M — 2.1 ± 0.2 M Q or of the mass-radius relation from the thermal 
emission of the isolated neutron star RX J1856-3754 pointing to either large radii of R > 14 km for a 
typical NS mass of < 1.4 M Q or large masses M > 2.0 M Q for radii not exceeding 12 km [12J. These 
measurements clearly demand a stiff equation of state and exclude standard models for hyperonic or 
quark matter interiors as well as mesonic condensates, see |13j . 

As has been argued in [T4], there are several modern QCD- motivated quark matter EsoS which 
could provide enough stiffness of high-density matter to be not in conflict with the new mass and 
mass-radius constraints. This leads, however, to the effect that hybrid stars with quark matter interiors 
"masquerade" as neutron stars [15] since they cannot be distinguished from each other by their mechan- 
ical properties, see also the contribution by T. Klahn in this volume [16] for a color superconducting 
three-flavor quark model with selfconsistently determined quark masses and pairing gaps [5*| [T7]. 

Therefore, also suggested signals from the timing behavior of pulsar spin-down [18J, frequency clus- 
tering [19] or population clustering [201 121] of accreting NSs should not be applicable. 

In the present work, we suggest a sensitive tool for "unmasking" the composition of neutron stars 
which is based on their cooling behavior. As the cooling regulators such as neutrino emissivities, heat 
conductivity and specific heat in quark matter might be qualitatively different from those in nuclear 
matter, due to the chiral transition and color superconductivity with some possibly sensible density 
dependence, the TA curves for hybrid stars could be significantly different from those of neutron stars. 
In order to reach the goal of unmasking the neutron star interior we introduce here a new method for the 
quantitative analysis of the cooling behavior consisting in the extraction of a NS mass distribution from 
the (yet sparse) TA data and its comparison with the (most likely) mass distribution from population 
synthesis models of NS evolution in the galaxy [22] . 

2 Cooling curves in the TA diagram 



The cooling behavior of compact stars belongs to the most complex phenomena in astrophysics. 
Therefore, the codes for its numerical simulation as developed by a few groups contain inputs (cooling 
regulators) of rather different kind, see, e.g., [231 [2U (2S1 12E1 [21] • Attempts to develop a Minimal Cooling 
Paradigm [21] by omitting important medium effects on cooling regulators [26| [28] unfortunately result 
in inconsistencies and suffer therefore from the danger of being not reliable. To develop a paradigmatic 
cooling code as an open standard, however, is rather necessary to cross-check the present knowledge 
of the groups before more sophisticated mechanisms like anisotropics due to the magnetic field [23] or 
special processes in the NS crust or at the surface are taken into account. Therefore, it is still premature 
to attempt an identification of the NS interior from the cooling behavior. 

In order to circumvent such a model dependence we employ a given cooling code developed in Refs. 
[251 [2S] and vary the matter properties such as EoS, superconductivity and star crust model such as to 
fulfill all constraints known up to now (mass, mass-radius, TA, brightness, etc.). Moreover, we try to 
use consistent inputs. 

We consider the cooling evolution of young neutron stars with ages t ~ 10 3 — 10 6 yr which is governed 
by the emission of neutrinos from the interior for t ^ 10 5 yr and thermal photon emission for t Z 10 5 yr. 
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Figure 1: Density dependence of the pairing gaps in nuclear matter together with that of the hypothetical 
X-gap in quark matter (left). Mass-central density relation for the two hadronic EoS models DBHF and DD-F4 
(right). The dot indicates the onset of the DU process. 



The internal temperature is of the order of T ~ 1 keV. This is much smaller than the neutrino opacity 
temperature T opac ~ 1 MeV as well as critical temperatures for superconductivity in nuclear (T c ~ 1 
MeV) or quark matter (T c ~ 1 — 100 MeV). Therefore, the neutrinos are not trapped and the matter 
is in a superconducting state. In Fig. [1] we show the density dependence of the pairing gaps in nuclear 
matter [301 122] together with that of the hypothetical X-gap in quark matter [311 [321 [33] . The phase 
transition occurs at the critical density n c = 2.75 n Q = 0.44 fm~ 3 . 

The main neutrino cooling processes in hadronic matter are the direct Urea (DU), the medium 
modified Urea (MMU) and the pair breaking and formation (PBF) whereas in quark matter the main 
processes are the quark direct Urea (QDU), quark modified Urea (QMU), quark bremsstrahlung (QB) 
and quark pair formation and breaking (QPFB) [34J. Also the electron bremsstrahlung (EB), and the 
massive gluon-photon decay (see [32]) are included. 

The lSo neutron and proton gaps in the hadronic shell are taken according to the calculations by [30] 
corresponding to the thick lines in Fig. 5 of Ref. [25J. However, the 3P2 gap is suppressed by a factor 
10 compared to the BCS model calculation of [30], consistent with arguments from a renormalization 
group treatment of nuclear pairing [36J . Without such a suppression of the 3P2 gap the hadronic cooling 
scenario would not fulfill the TA constraint, see [28J. 

The possibilities of pion condensation and of other so called exotic processes are included in the 
calculations for purely hadronic stars but do not occur in the hybrid ones since the critical density for 
pion condensation exceeds that for deconfinement in our case [25] . While the hadronic DU process occurs 
in the DBHF model EoS for all neutron stars with masses above 1.27 M & , it is not present at all in the 
DD-F4 model, see the right panel of Fig. [TJ We account for the specific heat and the heat conductivity 
of all existing particle species contributing with fractions determined by the f3— equilibrium conditions. 
Additionally, in quark matter the massless and massive gluon-photon modes also contribute. 

In the 2SC phase only the contributions of quarks forming Cooper pairs (say red and green) are 
suppressed via huge diquark gaps, while those of the remaining unpaired blue color lead to a so fast 
cooling that the hybrid cooling scenario becomes unfavorable [32]. Therefore, we assume the existence 
of a weak pairing channel such that in the dispersion relation of hitherto unpaired blue quarks a small 
residual gap can appear. We call this gap Ax and show that for a successful description of the cooling 
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scenario Ax has to have a density dependence. We have studied the ansatz 



A x = A exp 



-a 



A* - He 

, He 



where \i is the quark chemical potential, \i c = 330 MeV. For the analyses of possible models we vary 
the values of a and A , given in the Table 1 of [33] and shown in the left panel of Fig. [TJ 

The physical origin of the X-gap is not yet identified. It could occur, e.g., due to quantum fluctuations 
of color neutral quark sextet complexes [3]. Such calculations have not yet been performed within the 
relativistic chiral quark models. The size of the small pairing gaps in possible residual single color/ 
single flavor channels [37] is typically in the interval 10 keV - 1 MeV, see discussion in [38] . The specific 
example of the CSL phase is analyzed in more in detail in Refs. j39j HO, I4T] . 




Figure 2: Hadronic star cooling curves for DBHF model EoS. Different lines correspond to compact star mass 
values indicated in the legend (in units of M Q ), data points with error bars are taken from Ref. [24] . 



In Fig. 2 we present TA-diagrams for two hadronic star EsoS, where DBHF is an ab-initio calcula- 
tion for the Bonn-A nucleon-nucleon potential within the Dirac-Brueckner-Hartree-Fock approach [42J, 
discussed in the context of compact star constraints in [43, 44J. DD-F4 denotes a relativistic mean- 
field model of the EoS with density- dependent masses and coupling constants adjusted to mimick the 
behavior of the DBHF approach H6] . In the hadronic cooling calculations presented in Fig. 2 the 
crust model has been chosen according to Tsuruta's formula for the T m — T s relationship between be 
temperatures of the inner crust and the surface [23] . 

In Fig. 3 we show the TA diagrams for two hybrid star cooling models presented in Ref. [33J. 
The TA data points are taken from [24J. The hatched trapeze-like region represents the brightness 
constraint (BC) [17J. For each model nine cooling curves are shown for configurations with mass values 
corresponding to the binning of the population synthesis calculations explained in [33J . 

For the hybrid cooling scenario in [33J a more detailed measure for the ability of a cooling model to 
describe observational data in the temperature-age diagram had been introduced. Also the logN-LogS 
distribution constraint has been considered. In the TA diagram to encode the likelihood that stars in 
that mass interval can be found in the solar neighborhood, in accordence with the population synthesis 



4 



scenario, see Fig. [3] we use a marking with five grey values. The darkest grey value, for example, 
corresponds to the most populated mass interval 1.35 - 1.45 M & predicted by the mass spectrum used 
in population synthesis. 




4 5 6 7 2 3 4 5 

log 10 (t[yr]) log 10 (t[yr]) 



Figure 3: Cooling curves for hybrid star configurations with 2SC+X pairing pattern and X-gap model I (left) 
versus model IV (right). For the gaps see the left panel of Fig. 1. The grey value for the shading of the mass 
bin areas corresponds to the probability for that mass bin value in the population synthesis model of Ref. [22J. 



3 Mass distribution from TA data 

In the present work we suggest a method to construct a NS mass distribution as a tool for the quantita- 
tive characterisation of a given cooling model using the yet sparsely distributed measurements of cooling 
neutron stars in the TA plane as a probability measure. This method is described as follows. Once the 
cooling model is defined by the EoS and the cooling regulators, the configuration corresponding to a 
gravitational mass M and its cooling curve T(t; M) can be determined. Next a set of mass values Mj, 
i — 0, . . . , Nm is introduced, defining the borders of Nm mass bins as, for example, in the population 
synthesis. A pair of neighboring cooling curves T(t; Mj) and T(t; Mj_i) corresponds to the i th mass bin 
and delimits a strip in the TA plane. We now construct a measure for the number of cooling objects to 
be expected within this mass bin 

^Vcooi . /.T(t;Mi_i) 

Ni= J2 dt dTPj(T,t), (2) 

|— [ J JT{t;Mi) 

where iV cool denotes the total number of observed coolers used for the analysis and Pj (T, t) is the 
probability density to find the j th object at the point (T, t) in the TA plane. For the present exploratory 
study we make the simplest ansatz that Pj(T, t) is constant in the rectangular region defined by the upper 
and lower limits of the confidence intervals corresponding to the temperature and age measurements, 
(Tji,T ju ) and (tji,t ju ), respectively, 

P 3 (T,t) = [{t jt - t ju ){T jt - 7} U )] _1 6(T - T 3l )Q{T 3U - T)Q{t - t 3l )Q{t JU - t) . (3) 
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Note that in the case when the exact age the object is known (e.g., for a historical supernova), the 
time-dependence of Pj(T,t) degenerates to a 5-function and the t-integral in (TS]) can be immediately 
carried out, leaving us with a one-dimensional probability measure. 

This method has been applied to the cooling models for hadronic and hybrid stars described in 
the previous section. The results for the extracted mass distributions are normalized to 100 objects, 
defining N(M) = 100 ^/(E^q Nt), and shown in Fig. 4. 

As we see from Fig. 4, the results are very sensitive to the chosen cooling model. In the hadronic 
scenario the onset of the DU cooling mechanism drastically narrows the mass distribution around the 
critical mass for the DU onset, see Fig. 1. On the other hand the slow cooling model predicts more 
massive objects than could be justified from the independent population analysis. 

When comparing the density dependence of the pairing gaps, given in the left panel of Fig. 1, with 
the extracted mass distributions for the corresponding hybrid models in the right panel of Fig. 4, the 
direct relationship between the superconductivity and the mass distribution becomes obvious. 

The DU problem as it was previously discussed in the literature [251II31II8] was based on the intuitive 
understanding that the mass distribution can not be peaked at a critical mass value which accidentally 
is unique for all observed young objects. Our modification of the definition of the DU problem does not 
contradict that suggestion, but rather provides an additional measure which rehabilitates the validity 
of cooling scenarios including the DU process. 

On the other hand, the EoS model should obey the mass constraints too. Therefore, using the 
models discussed in this work we demonstrate that the most preferable structure of the compact object 
is likely to be a hybrid star with properly defined color superconductivity of the quark matter state in 
the core. 
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Figure 4: NS mass spectra extracted from the distribution of cooling data for both hadronic EoS models (left 
panel) and for hybrid stars with X-gap models I- IV (right panel). For comparison, the mass distribution of 
young, nearby NS from the population synthesis of Popov et al. [33] is shown at the bottom of the panels. 
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4 Conclusion 



In this contribution we have presented mass distributions obtained by analysing cooling calculations for 
NSs and demonstrated that they provide a sensible measure for the composition of compact star matter, 
even if the mechanical properties of the compact objects are almost identical. The comparison with 
a mass distribution from a population synthesis approach allows to favor hybrid stars with properly 
defined color superconducting quark matter core over other hybrid star or pure neutron star models. 
The presented approach makes testable predictions and will be quantitatively improved when a better 
statistics for TA data will emerge from present and future observational programmes. On the other 
hand, an improvement of the population synthesis might be required in connection with the ongoing 
development of supernova modeling and a deeper understanding of astrophysical processes in the galactic 
neighborhood [49] . The new method presented here to extract a NS mass distribution from the cooling 
behavior has proven useful as a new tool to unmask the NS interior since it is sensitive to subtle changes 
in the cooling regulators. It shall therefore be further developed and not be missed out when Astronomy 
meets QCD [50J and more stringent tests are applied to the behavior of nuclear matter under conditions 
of extreme densities in neutron stars, a place where the nature of the deconfinement transition can be 
studied. 
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